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Abstract. This paper has two agendas. Firstly, we exhibit new results for cov- 
erage processes. Let p(n, m, a) be the probability that n spherical caps of angular 
radius a in S m do not cover the whole sphere S m . We give an exact formula for 
p(n, m, a) in the case a € [tt/2, tt] and an upper bound for p(n, m, a) in the case 
a £ [0, 7r/2], which tends to p(n,m, 7r/2) when a — » n/2. In the case a £ [0, 7r/2] 
this yields upper bounds for the expected number of spherical caps of radius a 
that are needed to cover 5"". 

Secondly, we study the condition number ^(A) of the linear programming feasi- 
bility problem 3x G R m+1 Ax < 0, x ± where A £ R" x ( m + 1 ) i s randomly chosen 
according to the standard normal distribution. We exactly determine the distri- 
bution of ^(A) conditioned to A being feasible and provide an upper bound on 
the distribution function in the infeasible case. Using these results, we show that 
E(ln'^'( J 4)) < 21n(m + 1) + 3.31 for all n > m, the sharpest bound for this ex- 
pectancy as of today. Both agendas are related through a result which translates 
between coverage and condition. 
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1 Introduction 



1.1 Coverage processes on spheres 

One of the oldest problems in the theory of coverage pro- 
cesses is that of calculating the chance that a given region is 
completely covered by a sequence of random sets. Unfortu- 
nately there is only a small number of useful circumstances 
where this probability may be calculated explicitly. 

P. Hall [13, §1.11] 

In 1897 Whitworth [30] considered the following problem. Assume we place n 
arcs of angular radius a in the unit circle S 1 , whose centers are independently and 
randomly chosen from the uniform distribution in S 1 . What is the probability 
that these arcs do not cover S 1 ? 

Whitworth's problem is arguably at the origin of the theory of coverage pro- 
cesses. It was not until 1939 that an answer to the problem was given, when 
Stevens [28] showed that the probability in question is 

j =1 

where k = [^\ . Extensions of this result to other quantities related with random 
arcs in S* 1 are given in [25] . Extensions to random arcs with different lengths are 
given in [9, 16] and in [26] where an exact formula for the probability above is 
given for randomly placed arcs having random independent size. 

The extension of the original problem in S 1 to the 2-dimensional unit sphere 
S* 2 was considered by Moran and Fazekas de St. Groth [20]. Let p(n, a) denote 
the probability that n spherical caps of angular radius a and centers randomly 
and independently chosen from the uniform distribution on S 2 do not cover S 2 . 
Moran and Fazekas de St. Groth exhibited an approximation of p(n,a), and 
numerically estimated this quantity for a = 53° 26' (a value arising in a biological 
problem motivating their research). Shortly thereafter, Gilbert [10] showed the 
bounds 

(1 - A)" < p(n, a)<^ n(n - 1)A(1 - A)"" 1 , (2) 

where A = (sin^) 2 = ^(1 — cos a) is the fraction of the surface of the sphere 
covered by each cap. In addition, Gilbert conjectured that, for n — > oo, p(n,a) 
satisfies the asymptotic equivalence 

p(n, a) wn(n- 1)A 2 (l - A 2 )™" 1 . 
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This conjecture was proven by Miles [18], who also found an explicit expression 
(cf. [17]) for p(n,a) if a € [ir/2,ir], namely 



p{n, a) = Q jT a s in 2 ( n - 2 )(0/2)sin(2^)^ 

+ ^ Q jT a sm 2( - n -V (9/2) sin 3 6d9. (3) 

More on the coverage problem for S* 1 and S 2 can be found in [27]. Extensions of 
these results to the unit sphere S m in M m+1 for m > 2 are scarse. Let p(n, m, a) 
be the probability that n spherical caps of angular radius a in S m do not cover 
S m . That is, for a £ [0, 7r], and a±, . . . , a n randomly and independently chosen 
points in S m from the uniform distribution, define 

p(n, m, a) := Prob {S m / U™ =1 cap(aj, a)} , 

where cap(a, a) denotes the spherical cap of angular radius a around a. It can 
easily be seen that for n < m + 1 and a < tt/2 we have p(n, m, a) = 1. Moreover, 
Wendel [29] has shown that 

p(n,m,vr/2) = 2 1 ^f;( n 7 1 ). (4) 



k=0 



Furthermore, a result by Janson [15] gives an asymptotic estimate of p(n, m, a) 
for a — > 0. Actually, Janson's article covers a situation much more general than 
fixed radius caps on a sphere and it was preceded by a paper by Hall [12] where 
bounds for the coverage probability were shown for the case of random spheres 
on a torus. 

A goal of this paper is to extend some of the known results for S 1 and S 2 to 
higher dimensions. To describe our results we first introduce some notation. We 
denote by 

m + l 

O m :=vo\ m (S m ) 2T> 2 



r(afi) 

the m-dimensional volume of the sphere S m . Also, for t G [0,1], denote the 
relative volume of a cap of radius arccost e [0,7r/2] in S m by X m (t). It is well 
known that 

/>arccos t 

A m (t) = ^ zi / (smer-Ue. (5) 

v m Jo 

Our results are formulated in terms of a family of numbers C(m, k) defined for 
1 < k < m. These numbers are defined in §4.1 and studied in §5, where we give 
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bounds on C(m, k) and derive a closed form for k G {1, m — 1, m}. Furthermore, 
we will show that, for each m, the C(m, k) can be obtained as the solution of a 
system of linear equations, which easily allows to produce a table for their values 
(cf. Table 1). 

A main result in this paper is the following. 

Theorem 1.1 Let n > m > 1, a G [0, ir], and e = cos(7r — a). For a G [f ,vr] 

m { \ l' 1 
p(n,m,a)=J2[ k + 1 )C(m,k) J t m ~ k (1 - t 2 )^ 1 X^ 1 (t) dt 

and for a G [0, |) we have 

p(n, m, a) < T ^SP + ( ro + i) C ( m < ™) / ^ " l ^ i 1 ^W)""™" 1 *■ 

We remark that this formula, for a G [7r/2,7r] and m = 2, is identical with 
the one given by Miles (3). Also, for a G [0,7r/2] and m = 1, our upper bound 
for p(n, 1, a) coincides with the first term in Steven's formula (1) (cf. Remark 4.9 
below) . 

We may use Theorem 1.1, together with the bound on the C(m, k), to derive 
bounds for the expected value of N(m, a), the number of random caps of radius 
a needed to cover S m . The asymptotic behavior of N(m, a) for a — > has been 
studied by Janson [15]. Otherwise, we have not found any bound for E(N(m, a)) 
in the literature. 

Theorem 1.2 For a G (0, § ] we have 

E(N(m, a)) < 3m + 2 + y/m(m + 1) cos(a)A m (cos(a)) -2 ( - - 1 - 

V2A m (cos(a 

1.2 Polyhedral conic systems and their condition 

Among the number of interrelated problems collectively known as linear pro- 
gramming, we consider the following two. 

Feasibility of polyhedral conic systems (FPCS). Given a matrix A G 
W ix ^ m+1 \ decide whether there exists a nonzero x G R m+1 such that Ax < 
(componentwise) . 

Computation of points in polyhedral cones (CPPC). Given a matrix 
A g ]R nx ( m+1 ) such that y = {x G K m+1 | Ax < 0} / 0, find x G & . 
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By scaling we may assume without loss of generality that the rows ai,...,a n 
of A have Euclidean norm one and interpret the matrix A as a point in (S m ) n . 
We say that the elements of the set 

Fn,m ■= {A G (S m ) n \3xeS m <oi, x) < 0, . . . , (a n ,x) < 0} (6) 

are feasible. Similarly, we say that the elements in 

K,m ■= ^ € (S m ) n \3x £ S m (ai,x) < 0, . . . , (a n ,x) < 0} (7) 

are strictly feasible. Elements in (S m ) n \ T n>m are called infeasible. Finally, we 
call ill-posed the elements in S„ im := T n , m \ T° m . 

For several iterative algorithms solving the two problems above, it has been 
observed that the number of iterations required by an instance A increases with 
the quantity 

= \ 

dist(A, E n>m ) 

(here dist is the distance with respect to an appropriate metric; for the precise 
definition we refer to §2.1). This quantity, known as the GCC- condition number 
of A [11, 3], occurs together with the dimensions n and m in the theoretical 
analysis (for both complexity and accuracy) of the algorithms mentioned above. 
For example, a primal-dual interior-point method is used in [6] to solve (FPCS) 
within 

0(y/m + n(ln(m + n) + \n^(A))) (8) 

iterations. The Agmon-Motzkin-Schonberg relaxation method 1 [1, 21] or the per- 
ceptron method [23] solve (CPPC) in a number of iterations of order 0{ c t9{A) 2 ) 
(see Appendix B of [5] for a brief description of this). 

The complexity bounds above, however, are of limited use since, unlike n and 
m, 'tf(A) can not be directly read from A. A way to remove 'if (A) from these 
bounds consists in trading worst-case by average-case analysis. To this end, one 
endows the space ( 1 S" 11 ) ra of matrices A with a probability measure and studies 
^(A) as a random variable with the induced distribution. In most of these works, 
this measure is the unform one in (5"")™ (that is, matrices A are assumed to have 
its n raws independently drawn from the uniform distribution in S m ). 

Once a measure has been set on the space of matrices (and in what follows 
we will assume the uniform measure in (S m ) n ) an estimate on E(ln^(A)) yields 
bounds on the average complexity for (FPCS) directly from (8). For (CPPC) 

This method gives also the context in which 'if (A) was first studied, although in the feasible 
case only [11]. 
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the situation is different since it is known [5, Cor. 9.4] that E( < ^(^4) 2 ) = oo. Yet, 
an estimate for e > on 



Prob{^(^) >l/e\Ae T n . m ] 

yields bounds on the probability that the relaxation or perceptron algorithms 
will need more than a given number of iterations. Efforts have therefore been 
devoted to compute the expected value (or the distribution function) of ^(A) for 
random matrices A. 

Existing results for these efforts are easily summarized. A bound for E(ln 'if (A)) 
of the form C(min{n,mlnn}) was shown in [4]. This bound was improved [7] 
to maxjlnm, In Inn} + 0(1) assuming that n is moderately larger than m. Still, 
in [5], the asymptotic behavior of both ^(A) and In ^(A) was exhaustively stud- 
ied and these results were extended in [14] to matrices A £ (^g m ^ n drawn from 
distributions more general than the uniform. Independently of this stream of 
results, in [8], a smoothed analysis for a related condition number is performed 
from which it follows that E(ln^(A)) = O(lnn). 

Our second set of results adds to the line of research above. Firstly, we provide 
the exact distribution of ^(A) conditioned to A being feasible, and a bound on 
this distribution for the infeasible case. 

Theorem 1.3 Let A £ (g" 1 )™ be randomly chosen from the uniform distribution 
in (S m ) n , n> m. Then, for e G (0, 1], we have 



Pmbl 1g{A) > 1/e 



- X 

n,m f 

on-1 m / \ re 

Lfe=o I k ) fe =i V s + V Jo 



Prob\tf(A) > 1/e 



>n — 1 



< 



jk—m-\-l 



(V) v™ + 



n 



\C(m,m) (l-t 2 )—(l-\ m {t)) 



n—m—1 



dt. 



Secondly, we prove an upper bound on E(ln^(^4)) that depends only on m, 
in sharp contrast with all the previous bounds for this expected value. 



Theorem 1.4 For matrices A randomly chosen from the uniform distribution 
in (S m ) n with n>mwe have E(lntf(A)) < 21n(m + 1) + 3.31. 
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Note that the best previously established upper bound for F,(ln < ^'(A)) (for 
arbitrary values of n and m) was O(lnn). The bound 21n(m + 1) + 3.31 is not 
only sharper (in that it is independent of n) but also more precise (in that it does 
not rely on the O notation) 2 . 

1.3 Coverage processes versus condition numbers 

Theorems 1.1 and 1.3 are not unrelated. Our next result, which will be the first 
one we will prove, shows a precise link between coverage processes and condition 
for polyhedral conic systems. 

Proposition 1.5 Let a±,...,a n be randomly chosen from the uniform distri- 
bution in S m . Denote by A the matrix with rows a±, . . . ,a n . Then, setting 
e := | cos(a)| for a £ [0, n], we have 



While Proposition 1.5 provides a dictionary between the coverage problem 
in the sphere and the condition of polyhedral conic systems it should be noted 
that, traditionally, these problems have not been dealt with together. Interest on 
the second focused on the case of ^(A) being large or, equivalently, on a being 
close to 7r/2. In contrast, research on the first mostly focused on asymptotics for 
either small a or large n (an exception being [29]). 

Acknowledgments. We thank Dennis Amelunxen for pointing out an error in 
a previous version in the paper. 

2 Main Ideas 

In this section we describe in broad strokes how the results presented in the 
introduction are arrived at. In a first step (§2.1), we give a characterization of 
the GCC condition number which establishes a link to covering problems, thus 
leading to a proof of Proposition 1.5. We then proceed by explaining the main 
ideas behind the proof of Theorem 1.3. 

In all that follows, we will write [n] = {1, . . . , n} for n G N. 

2 Recently a different derivation of a O(lnm) bound for E(ln'rf(A)) was given in [2]. However, 
this derivation does not provide explicit estimates for the constant in the O notation. 




yLA + p ro b{A(jLT n , m andtf{A) > ±} if a e [0,tt/2]. 
/2) = Prob{A G F n , m } = 2 1 - IXo (V) • 




In particular, p(n, m, n 
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Figure 1: A SIC with a e (0,tt/2) (left) and with a e (n/2,n) (right). 



2.1 The GCC condition number and spherical caps 

A key ingredient in what follows is a way of characterizing the GCC condition 
number in terms of spherical caps. For p G S m and a G [0, 7r], recall that 

cap(p, a) := {x G S 1 " 1 | (p, x) > cos a}. 

A smallest including cap (SIC) for 4 = (ai, . . . , a n ) G (5 m ) ra is a spherical cap of 
minimal radius containing the points ai, . . . ,a n . If p denotes its center, then its 
blocking set is defined as {i G [n] \ (p, a«) = cos a}, which can be seen as a set of 
"active" rows (cf. Figure 1). 

A largest excluding cap (LEC) for A is the complement of a smallest including 
cap for A. Note that (by a compactness argument) a SIC always exists, and while 
there may be several SIC for A, its radius is uniquely determined. For the rest 
of this article, we denote this radius by p(A) and set t(A) := cos p(A). The 
following is one of many equivalent ways [3, 5] of defining the GCC condition 
number. 

Definition 2.1 The GCC condition number of A G (S m ) n is defined as tf(A) := 
1/| cos ,o(j1)| G (l,oo]. 

In order to understand the relation of this condition number to distance to 
ill-posedness we review a few known facts (for more information see [3] and [5]). 
Recall the definition of T n , m and T^ m given in Equations (6) and (7). It is easy 
to see that J- n , m is a compact semialgebraic set with nonempty interior T^m- 
The set S njm := J- nm \ J~nm ls the- topological boundary of J- n ,m- It consists of 
the feasible instances that are not strictly feasible. Note that if n > m + 1, then 
S„ im is also the boundary of the set of infeasible instances T n ,m '■= {S m ) n \J r n , m . 
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Hence in this case S„ im consists of those instances that can be made both feasible 
and infeasible by arbitrarily small perturbations. 

The next lemma summarizes results from [3, Thm. 1] and [5, Prop. 4.1]. It 
is enough to prove Proposition 1.5. 

Lemma 2.2 Wehavep(A) < it /2 if and only if A £ T° im . Moreover, p(A) = ir/2 
if and only if A £ S njjn . □ 

Proof of Proposition 1.5. We claim that 

p(n, m, a) = Prob{p(A) < it — a}. (9) 

Indeed, U^ =1 cap(aj, a) ^ S m iff there exists y £ S m such that y cap(aj,a) for 
all i. This is equivalent to 3y\/iai cap(y,a), which means that an LEC for A 
has angular radius at least a. This in turn is equivalent to p{A) < ir — a, thus 
proving the claim. 

Equation (9) for a = ir/2 combined with Lemma 2.2 and Wendel's result [29] 
stated in Equation (4) yields 




= p{n,m,ir/2) = Prob{p(A) < vr/2} = Prob{A £ F n , m }. 



Suppose now a £ [ir/2,ir]. Then 

p(A) <ir-a p(A) < vr/2 and tf(A) < ^, 

showing the first assertion of Proposition 1.5. Furthermore, for a £ [0, |] 
p(A) < 7T - a <=^ p(A) < tt/2 or {^(A) > ir/2 and | cos p(A)\ < | cos(7r - a)|) , 
showing the second assertion of Proposition 1.5. □ 



2.2 Towards the proof of Theorem 1.3 

To prove the feasible case in Theorem 1.3 we note that 



1 



Prob <^ <jf (A) > - 



A £ J- n < 



1 



MO\T n . 



-vo\T n , m {e), 



where T n ^ m {e) = {A £ T° n m \ t(A) < e}. But vo\J : njm is known by Proposi- 
tion 1.5. Therefore, our task is reduced to computing vo\J 7 1hm (e). As we will 
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see in Section 3.1, the smallest including cap SIC(A) is uniquely determined for 
all A £ J-° m . Furthermore, for such A, t(A) depends only on the blocking set 
of A. Restricting to a suitable open dense subset 1Z n ,m{£) ^ ^n,m{ £ ) of "regular" 
instances, these blocking sets are of cardinality at most m + 1. This induces a 
partition 

n nim (e) = \Jn I nim (e), 

I 

where the union is over all subsets / C [n] of cardinality at most m + 1, and 
^-n,m( £ ) denotes the set of matrices in TZ n ^ m (e) with blocking set indexed by /. By 
symmetry, volTZ^ m (e) only depends on the cardinality of /, hence it is enough to 
focus on computing vol TZn^ {s) for k = 1, . . . , m. The orthogonal invariance and 
the particular structure of the TZn^ (s) (involving certain convexity conditions) 
makes possible a change of coordinates that allows to split the occurring integral 
into an integral over t and a quantity C(m, k) that depends only on m and k: 

vol TZ^ (e) = C(m, k) f g(t, n, m, k)dt. 

Jo 

More precisely, we proceed as follows: 

(1) By Fubini, we split the integral over TZn^\e) into an integral over the first 
k + 1 vectors a\, . . . , au+i (determining the blocking set [k + 1]) and an integral 
over afc + 2, • • • , a n taken from SIC^): 

voi^+%)=/ (f dis-r-"- 1 ) dn^i 

= f [k+1]/ G(A)dU^ m . (10) 

This is an integral of the function G(A) := vol (cap(p(A), p(A))) n ^ k ^ 1 , which is 
a certain power of the volume of the spherical cap SIC(vl). 

(2) The next idea is to specify A = (ai, . . . , a k +i) in Tl k k ^ m {e) by first specifying 
the subspace L spanned by these vectors and then the position of the on the 
sphere S m n L = S k (cf. Figure 2). 

Let Gfc + i(M m+1 ) denote the Grassmann manifold of (k + l)-dimensional sub- 
spaces in R m+1 and consider map 

^KL( £ ) - G k+1 (R m+1 ), (ai, . . . ,o fc+ i) ^ L = span{ai, . . . ,a k+1 }. (11) 

Clearly, a vector a 6 S m lies in the special subspace L$ := x iff it lies 

in the subsphere S k . Hence the fibre over Lq consists of all "regular" tuples 
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{S m ) k+1 G k+1 (W n+1 ) (S m n L = s k ) k+1 

Figure 2: Determining (01,02) € (S 2 ) 2 by first giving its span L e G2(K 3 ) and then 

a,eS 2 n L. 

A = (ai, . . . , cifc + i) G (iS*) fc+1 such that t(^4) < e, and hence the fibre can be 
identified with T^k+i\( £ )' Using the orthogonal invariance and the coarea for- 
mula (also called Fubini's theorem for Riemannian manifolds) we can reduce the 
computation of the integral (10) to an integral over the special fibre 
This leads to 

L +1] ( , G ^ d <"l = vd J G(A)J(A) dU^ k 

where J (A) is the normal Jacobian of the transformation (11). 

(3) To specify a regular A = (a±, . . . , a^+i) G (S k ) k+1 , we first specify the direc- 
tion p = p{A) G S k and the height t = t(A) G (0, e), and then the position of the 
di on the subsphere {a G S k \ (a,p) = t} ~ S ,fe ~ 1 (cf. Figure 3). 




Figure 3: Determining (ai,a 2 ,a 3 ) e (S 2 ) 3 by specifying the direction p, the height t, 
and the a, on the subsphere {a e S 2 = t} by b{. 
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More precisely, we consider the map 

- S k x (0 >e ), A~(p(A),t(A)). (12) 

The fibre over (p, t) , where po = (0, . . . , 0, 1) is the "north pole" , consists of tuples 
(ai . . . , cifc+i) lying on the "parallel" subsphere {a € S k \ (a,p) = t}. The vectors 
di can be described by points bi € S* -1 , which are obtained by projecting ai 
orthogonally onto R k x and scaling the resulting vector to length one. 

The orthogonal invariance and the coarea formula allow us to reduce the 
computation of the integral over TZ^^ k (e) to the integration over t £ [0, e] of an 
integral over the the special fibres over (po,t). The latter integral is captured by 
the coefficient C(m, k), which can be interpreted as a higher moment of the vol- 
ume of the simplex A spanned by random points bi, . . . , bu+i on the sphere S* -1 . 
However, we have to respect the convexity condition that the origin is contained 
in the simplex A spanned by the 6j, which complicated matters. Altogether, we 
are lead to a formula for volT?]^™ 1 ' (e) of the shape 

vo\G k+1 (R m+1 ) [ G(A)J(A)dR%£{ = C(m,k) f g(t,n,m,k)dt, 

where g(t, n, m, k) is obtained by isolating the part of the resulting integrand 
that depends on t. 

In order to implement this plan, we have to isolate the appropriate regularity 
conditions, i.e., to identify the sets TZ^ m (e), and to compute the normal Jacobians 
of the maps (11) and (12). For the latter task, we prefer to use the language of 
differential forms, as is common in integral geometry [24]. 

Unfortunately, the above argument above does not carry over to the infeasible 
case. Nevertheless, the ideas described above are sufficient to obtain the upper 
bound in Theorem 1.3. 

The rest of the paper proceeds as follows. In Section 3 we describe the ba- 
sic facts on smallest including caps and integration on manifolds that will be 
needed to make formal the ideas exposed above. Then, in Section 4, we prove 
Theorem 1.3. Theorem 1.1 immediately follows via Proposition 1.5. Finally, in 
Section 5, we give bounds for all, explicit expressions for some, and a way to com- 
pute the coefficients C(m,k). ^,From these bounds we derive the Theorems 1.2 
and 1.4. 
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3 Preliminaries 



3.1 Properties of smallest including caps 

Recall from Section 2.1 the definition of smallest including caps (SICs) for a 
given A = (ai, . . . , a n ) G (S m ) n . A crucial feature of our proofs is the fact that 
a strictly feasible A has a uniquely determined SIC. This is a consequence of the 
following crucial lemma, which provides an explicit criterion for a spherical cap 
being a smallest including cap of A. This lemma is a generalization of Lemma 4.5 



Lemma 3.1 (a) For a strictly feasible A G .F° m there exists exactly one small- 
est including cap. 

(b) Let (p,t) G S m x (0,1] and 1 < k < n. Suppose that (ai,p) = t for all 
i G [k + 1] and (o,i,p) > t for all i G [n] \ [k + 1]. Then cap(p, arccos t) is the 
smallest including cap of A if and only if 



Proof. We first show that Assertion (b) implies Assertion (a). Suppose cap(j>i, p) 
and cap(j>2;P) are SICs for A and put t := cos p. Note that £ > 0. Assertion (b) 
implies that tpi is contained in the convex hull of a±, . . . , a n , hence there exist Aj > 
such that ^ \ = 1 and tpi = \a-i- Therefore, (tpi,p2) = J2i ^i{ a uP2) > t, 
as (aj,p2) > t for all i. This implies (pi,P2) > 1 and hence p\ = p2- 

The proof of Assertion (b) goes along the lines of Lemma 4.5 in [5]. Suppose 
first that cap(p, a) is an SIC for A, where a := arccost. It is sufficient to 
show that p £ conejdi, . . . , a^+i}. Indeed, if p = Yl ^i a i with A« > 0, then 
tP = XX^iK- Furthermore, J2(t\i) = J2^i( a i,P) = (Y^i a i,p) = \\p\\ 2 = 1- 
Hence tp G convjai, . . . , a^+i}- 

We now argue by contradiction: if p is not contained in conejai, . . . , a^+i}, 
then there exists a vector v G S m such that (p,v) < and (a,i,v) > for 
alH G {1, . . . , k + 1}. For <5 > we set 



in [5]. 



G convjai, . . . 



p + Sv p + 5v 



(13) 



: = 



Then for 1 < i < fc + 1 and sufficiently small 5 we have 




where we used that (a,i,p) = t, (ai,v) > 0, and (p,v) < 0. 
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For k + 2 < i < n the function 5 — > (a,i,ps) is continuous at <5 = 0. Since, 
by hypothesis, (ai,p) = (ai,po) > t, it follows that {ai,p$) > t for 5 sufficiently 
small. From this we conclude that for sufficiently small 5 there exists t$ > t such 
that (ai,ps) > t$ for all i G [n]. Setting ay = arccost^ we obtain that as < a 
and aj G cap(p5,ay) for all i G [n], contradicting the assumption that cap(p, a) 
is a smallest including cap. 

To prove the other direction, we suppose tp G convjai, . . . , a^+i}. For q G S m 
let a(q) we denote the angular radius of the smallest spherical cap with center q 
containing ai, . . . a n . If we assume that cap(p, a) is not a SIC for A, then there 
exists a vector u G S m and <5o > 0, such that (v,p) = and, for all < 5 < So, 
a (ps) < a(p), where pg = (i-e., we have a direction v along which we can 

move to obtain a smaller cap). This means that 

min (oi,p s ) > min (a^p) = t. 

±<i<n l<i<n 

Therefore, for all i G [A; + 1] we have 

(ai,p) + «5(aj,w) . . 

<Oi,P«> = ^l + $2 > 1 = ( a " P < 

for sufficiently small 5, which implies that (a«, v) > 0. Let n G be such that 
tP = Ei<i<fe+i MiOi and Ei<j< fc +i Mi = 1- Then we have 

t(P, v ) = ^2 Vi(ai,v) > 0, 

l<i<fc+l 

contradicting the assumption that (p, v) = 0. Thus cap(p, a) is indeed a smallest 
including cap. □ 

For a strictly feasible A, we denote the center of its uniquely determined SIC 
by p(A) and its radius by p{A). The blocking set BS(A) of A is defined as the 
blocking set of the SIC of A. It is not hard to see that BS(^4) can have any 
cardinality greater than one. 

However, we note that in the infeasible case, there may be more than one 
smallest including cap: consider for instance three equilateral points on the circle 
(right-hand side in Figure 1). It is known [5, Prop. 4.2] that in this case, the 
blocking set of a SIC has at least m + 1 elements. In the infeasible case, one 
direction of the characterization of smallest including caps of Lemma 3.1 still 
holds. The proof is similar as for Lemma 3.1. 

Lemma 3.2 Let cap(p, arccos t) be a SIC for A G (S m ) n with p G S m and 
t G (—1, 0). Suppose (di,p) = t for i G [m + 1] and (di,p) > t for i = m + 2, . . . , n. 
Then tp G convjai, . . . , a m+ \}. 
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Proof. Suppose tp G" convjdi, . . . , a m+ i}. Then — p G" conejdi, . . . , a m +i} and 
hence there exists a vector v G S m such that {—p,v) < and (dj,t>) > for 
all i. For 5 > we define ps as in (13). Take 5 sufficiently small so that 
t < t + 5{a,i,v) < for all i 6 [m + 1]. Then, for i £ [m + 1] and 5 sufficiently 
small, we have 

t + 6(ai,v) 
(ai,ps) = . > t, 

^/l + 25(p,v} + 5 2 

where we used that (ai,p) = t, (ai,v) > 0, and (p,v) > 0. This shows that 
cap(p, arccos t) is not a smallest including cap. □ 

We present a few more auxiliary results that are needed for the proof of our 
main result. 

Lemma 3.3 For given linearly independent a±, ... ,0,^+1 £ S m , 1 < k < m, 
there exist unique p G S m and t G (0, 1) such that 

p G span{ai, . . . , afc+i} and (aj,p) = t for all i G [k + 1]. 

Moreover, the map (ai, . . . , 0^+1) 1— > (p, t) is differentiable. 

Proof. Let .4 denotes the affine span of a±, . . . , afc+i, L the underlying linear 
space, and C the linear span of A. Since the are linearly independent, we have 
A 7^ L and thus dim .A = dimL = k, dim£ = k + 1. Hence the intersection of 
C with the orthogonal complement L 1 - is one-dimensional and contains exactly 
two elements of length one. Take the one such that the common value t = (ai,p) 
is positive. This shows existence and at the same time the uniqueness of p, t. 

Suppose now k = m and let A denote the square matrix with the rows 
ai , . . . , d m+ i . The conditions (dj , p) = t can be written in matrix form as Ap = te, 
where e:=(l,...,l). By solving this equation we obtain the following explicit 
formulas 

p{A) = uk\\ A ~ le > t{A) = vh\\- (14) 

This shows the differentiability of the map A ^ (p, t) in the case k = m. We 
leave the proof in the general case to the reader. □ 

The next result, though very elementary, will be useful for clarification. 

Let p G S k and t / and consider elements a±,...,ak+i G S k satisfying 
(a,i,p} = t for all i. Let 6j G S' fc_1 be the scaled-to-one orthogonal projection 
of dj onto the orthogonal complement of M.p. That is, d« = rbi + tp, where 
r=(l-t 2 )V2. 

Lemma 3.4 The following conditions are equivalent: 
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1. The affine hull A of a±, . . . , at+\ has dimension k. 

2. The span of b±, ... , b^i has dimension k. 

3. oi, . . . , dk+i are linearly independent. 

Proof. The equivalence of the first two conditions is obvious. The equivalence 
of the first and third condition follows from dim spa n (A) = dim A + 1 (here we 
use t ^ 0). □ 



3.2 Volume forms on Grassmann manifolds 

Integration on Grassmann manifolds will play a crucial role in our proofs. We 
recall some facts about the relevant techniques from integral geometry and refer 
to Santalo's book [24, II. 9-12] and the article [19] for more information. We recall 
that volume elements are always unsigned forms. 

Let M be a Riemannian manifold of dimension m, p G M, and let y = 
(yi, . . . , y m ) T : U i ^ M m be local coordinates in a neighbourhood U of p such that 
d/dyi, . . . , d/dy m are an orthonormal basis of T p M. The natural volume form on 
M at p associated to its Riemannian metric is then given by dM = dy\f\ - ■ -Ady m , 
where dyi is the differential of the coordinate function yi at p. 

In the case of a sphere, we get such local coordinates around a point p G 
S m by projecting onto the orthogonal complement of p. More precisely, let 
(ei, . . . , e m +i) be an orthonormal basis of ]R m+1 satisfying e± = p (so that 
e2, • • • , e m +i span the tangent space T p S m ). For a point x = (x\, . . . , x m +i) T G 
S m in a neighbourhood of p set y« = (x, ei). Then (7/2, • • • , J/m+i) are local coor- 
dinates of S m around p such that d/dyi are pairwise orthogonal at p. Hence the 
volume element of S m at p is given by 

dS m = u> 2 A ■ ■ ■ Aai m+ i, 

where cjj := cZy, = (dx, e,) and cZx = (dxi, . . . , dx m +i) T . Hence, if we denote by 
E the (m+ 1) x (m+ l)-matrix having the as rows, we obtain the volume form 
by wedging the non-zero entries of Edx. 

In a similar fashion we define volume forms on Stiefel manifolds (for details 
and further justification we refer to [24]). A k- frame is a set of k linearly indepen- 
dent vectors. For 1 < k < m + 1, the Stiefel manifold Vfc(M m+1 ) is defined as the 
set of orthonormal /c-frames in M m+1 . It is a compact Riemannian submanifold 
of (S m ) k . Let Q = (qi, ...,%) G Vfc(M m+1 ) and (ei, . . . , e m+1 ) an orthonormal 
basis of IR m+1 , such that e\ = q\,...,ek = qt- Then the volume element of 



16 



V k (R m+1 ) at Q is given by 

dV k (R m+1 )= f\ K, + iA-A Wl , m+ i), 
i<i<k 

where Uij = (dqi, ej) for 1 < i < k and 1 < j < m+1. (In terms of the (m+1) x /c 
matrix dQ, this corresponds to wedging the entries below the main diagonal.) 
With this volume element we have volVfc(R m+1 ) = O m • • • O m+ i^k- 

We denote by Gfc(IR m+1 ) the Grassmann manifold of /c-dimensional subspace 
of IR m+1 . One way of characterizing it is as a quotient of a Stiefel manifolds, by 
identifying frames that span the same subspace. Let L 6 Gfc(M m+1 ) and choose 
a frame Q £ V/%(R m+1 ) spanning L. If V k (L) denotes the Stiefel manifold of 
orthonormal /c-frames in L and dVk(L) its volume element at Q, then it is known 
that the volume element dGk(R m+1 ) of the Grassmann manifold at L satisfies 
(see [19, (10)]) 

dV k (R m+1 ) = dG k (R m+1 ) A dV k (L). (15) 
^From this equality it follows that 

dG k (R m+1 )= f\ K W A-A Wj , m+ i), 
i<i<k 

with the Ldij as defined in the case of the Stiefel manifold. (In terms of the matrix 
E dQ, this corresponds to wedging the elements in the lower (m + 1 — k) x k 
rectangle.) As a consequence of (15), the volume of the Grassmannian is given 
by 

Gk,m+l '■= VOlGfc v - 

■ ■ ■ 

Equation (15) has a generalization to frames that are not orthogonal, i.e., to 
points in a product of spheres (S m ) k . Let L € G fc (M m+1 ) and set S(L) := LnS m , 
so that S(L) = S k ~ 1 . Choose a basis a\, . . . , of L consisting of unit length 
vectors, i.e., a point A = (ai, . . . , afc) in S(L) k . We denote by vo\(A) the volume 
of the parallelepiped spanned by the aj. Then the volume form of (S m ) k at A 
can be expressed as 

d(S m f = vo\(A) m - k+1 dG k (R m+1 ) A dS{L) k . (16) 

This equation can be derived as in [19], see also [24, II. 12.3]. It also follows 
as a special case of a general formula of Blaschke-Petkantschin type derived by 
Zahle [31], see also the discussion in [22]. 

A beautiful application of Equation (16) is that it allows an easy computation 
of the moments of the absolute values of random determinants. The following 
Lemma is an immediate consequence of (16) (see also [19]). 



pm+l\ _ @m+l~k "'Or, 
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Lemma 3.5 Let B £ (S k ) k+1 he a matrix with rows bi, . . . ,bk independently 
and uniformly distributed in S k . Then 

v(\det(B) r -^) = (-^Y -J— . 
\ / yk,m+i 

4 The probability distribution of ^(A) 

This section is devoted to the proof of Theorem 1.3. 
4.1 The feasible case 

Recall that, for A £ 3~n,m-> we denote the center and the angular radius of the 
unique smallest including cap of A by p(A) and p(A), respectively, and we write 
t(A) = cos p(A). 

Our goal here is to prove the first part of Theorem 1.3, for which, as we noted 
in §2.2, we just need to compute the volume of the following sets 

F n , m (e) := {A G F° n>m \ t(A) < e}. 

For this purpose it will be convenient to decompose J- nt m(e) according to the size 
of the blocking sets. Recall that the blocking set of A G T° i m is defined as 

BS(A) = {ie[n]\{p(A),a i )=t(A)}. (17) 

For / C [n\ with |/| < n and e G (0, 1] we define ^ m (e) to be the set of all 
A £ F n ,m{e) such that BS(A) = /. 

For technical reasons we have to require some regularity conditions for the 
elements of F^mi 6 )- 

Definition 4.1 We call a family (a±, . . . , a^i) of elements of a vector space 
centered with respect to a vector c in the affine hull A of the dj if dim^l = k 
and c lies in the relative interior of the convex hull of the Oj. We call the family 
centered if it is centered with respect to the origin. We now define, for / C [n], 

^■n,m( £ ) := £ Fn,m( £ ) I { a i)i&i ls centered with respect to j. 
Note that, by definition, the a« are aflinely independent if |/| < m + 1. 
Lemma 4.2 1. ^ m (e) is of measure zero if \I\ > m + 1. 
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2. If \I\ < m + 1, then IZ^ ^e) is open in (S m ) n and J r n,m{ £ ) JS contained in 
the closure ofTZ^^e). 

3. J 7 n,m( £ ) \ ^-n,m( £ ) nas measure zero. 

Proof. 1. If A G ^n,m( £ )^ then {aj | i G /} is contained in the boundary of the 
SIC of A and hence its affine hull has dimension at most m. On the other hand, 
the affine hull of (aj)j e / is almost surely R m+1 if |/| > m + 1. 

2. The fact that TZ^ m (e) is open in (S m ) n easily follows from the continuity 
of the map A i— > (p(A), t(A)) established in Lemma 3.3. 

Suppose now A G Fn,m{ £ )- By Lemma 3.1 we have t(A)p(A) £ convjaj | i G 
/} for A G J^ m (e). It is now easy to see that there are elements A' arbitrarily 
close to A such that A' is centered with respect t{A')p{A'). This shows the second 
assertion. 

3. By part two we have VJ n m {e) C ^[ m (e) C r R, I n m [e). Since we are dealing 
with semialgebraic sets, the boundary of TZ^ m (e) is of measure zero. □ 

It is clear that the T^ m wrt h / of the same cardinality just differ by a per- 
mutation of the vectors. Using Lemma 4.2 we obtain 



V0l^ n , m (£)= v O | ^n,m(e)=E( fc + l) VOI <m(e)- (18) 

|/|<m+l fc=l ^ ' 

where we have put lZ k m (e) := TZn^n\ £ ) to ease notation. 

Hence it is sufficient to compute the volume of TZ^ m (e). For this purpose we 
introduce now the coefficients C(m, k). 

Definition 4.3 We define for 1 < k < m 

(h\\m-k+l r 

C(m, k) = G k , m / (vol fc A) m - fc+1 d^*" 1 )^ 1 , 

where is the following open subset of S k ~ 1 

M k := {(6i,...,6 fc+ i) G (S k ~ 1 ) k+1 | (6i,...,6 fc+ i) is centered} 

and A : M k — > R maps I? = (&i, . . . , bk+i) to the convex hull of the fy. 

Example 4.4 We compute C(m, 1). Note that M^i = {(-1, 1), (1, -1)} and 
Q\,m = \O m -\- Hence 



C(m, 1) = -j- Jo m _! / (voliA) m dMi,i = 







m—1 
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The assertion in Theorem 1.3 about the feasible case follows immediately 
from Proposition 1.5, Equation (18), and the following result. 

Proposition 4.5 Let e G (0, 1]. The relative volume of Fnmi 6 ) JS given by 



C(m,k) [ t m - k (l - t 2 )^ km - 1 X m ( y t) n ' k - 1 dt. 
Jo 



Qn 

Proof. Consider the projection 

^n,m( £ ) ^fc+l,m( £ )> (°1 >•••>%) ^ ^ = («1 , • • ■ , Ofc+l)- 

By Lemma 3.1, this map is surjective and its fiber over A consists of all (A, a k +2, ■ ■ ■ , 
such that aj lies in the interior of the cap cap(p(A) , p(A)) for all « > k + 1. By 
Fubini, and using (5), we conclude that 



171 — k- 



T-= I X m (t(A)r- k - 1 d(S m ) k+1 . (19) 



We consider now the following map (which is well defined, cf. Lemma 3.4) 

ftfc+i,m( £ ) ~* G fe+ i(M m+1 ), (oi, . . . , Ofc+i) ^ L = span{ai, . . . , a k+1 }. 

We can thus integrate over T^+imi 8 ) by ^ rst integrating over L G Gfc+i(R m+1 ) 
and then over the fiber of L. By Equation (16), the volume form of (S m ) k+1 at 
A can be written as 

d(S m ) k+1 = vo\{A) m - k dG k+1 (R m+1 )(L) A dS(L) k+1 , 

where S(L) k+1 denotes (fc + l)-fold product of the unit sphere of L. By invariance 
under orthogonal transformations, the integral over the fiber does not depend on 
L. We may therefore assume that L = R fc+1 , in which case the fiber over L can 
be identified with TZ k +1 k (e). Thus we conclude from Equation (19) that 

VO '^ri £) = W / voKA)™-*^^^))"-*- 1 ^*)^ 1 . (20) 

In a next step, we will perform a change of variables in order to express the 
integral on the right-hand side of Equation (20) as an integral over t involving 
the coefficients C(m,k). 

Note that by Lemma 3.3, p(A) 6 S k and t(A) G (0, 1) are defined for any 
A £ GL(fc+l, R) and depend smoothly on A. A moment's thought (together with 
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Lemmas 3.1 and 3.4) shows that we have the following complete characterization 

(e) = {AG (S k ) k+1 | A is centered with respect to t(A)p(A), 
0<t(A)<e, Vi {oi,p(A))=t(A)y 

For A e K k k+lk (e) set P '■= P( A )> 1 '■= *(^)> and r := r (*) := C 1 ~ t 2 ) 1/2 - For 
i G [k + 1] we define bi as the scaled-to-one orthogonal projection of onto the 
orthogonal complement of Rp, briefly ai = rb{ + tp. The matrix B = B(A) with 
the rows b\, . . . , 6^+1 can be written as B = £ (A — tep T ). Clearly, B is centered. 
We define now 



W, 



k = G (.S^)^ 1 x 5 fc | Bp = and 5 is centered} . 

This is a Riemannian submanifold of (S k ) k+2 of dimension k(k + 1) — 1. We thus 
have a map 

<p k : n k k+hk {e) -^W k x (0,s), A^(B(a),p(A),t(A)). 

The inverse of this map is given by (B,p,t) i— > A = rB + tep T . It is well- 
defined since, by Lemma 3.4, A is invertible when B is centered. The Jacobian 
J (A) of the diffeomorphism (p k is stated in the next lemma, whose proof will be 
momentarily postponed. We remark that this lemma can also be derived from 
[22, Lemma 1] (a special case of Zahle's theorem [31]), with K being the unit 
ball. 

Lemma 4.6 The volume form of (S k ) k+1 at A can be expressed in terms of the 
volume form ofW k x (0, e) as follows 

h u., r (fc-2)(fc+i) yo \(A) 
d(S k ) k+1 = J(A) dW k A dt = — dW k A dt. 

We now express the Jacobian J (A) in terms of (B,p,t). The volume vol(A) 
of the parallelepiped spanned by the ctj equals {k + 1)! times the volume of the 
pyramid with apex and base A(ai, . . . ,a k+ ±), the latter denoting the simplex 
with vertices a\, . . . ,a k+ \. Moreover, this pyramid has height t and it is well 
known that the volume of a (k + l)-dimensional pyramid with height t and base 
B is j^j- times the (/c-dimensional) volume of B. This implies 

vo\(A) = (k+ 1)! ^ vol fc _iA(r6i, . . . ,rb k+1 ) = k\ r k tvo\ k A(B). 
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/,From this expression, together with (20), we conclude that 

wl 5 m(e) = gfc+ & +1 I voi(Ar-w rW(fc " 2)(fc+1) w*- 1 ^ fe ) fc+1 

— / vol fc A(B) m - fc+1 dVK fe / t m - fc r(t) fcm - 2 A m (t)"- fc - 1 dt. 



Consider the projection pr: W k — > S k ,(B,p) ^ p. We note that its fiber over 
p can be identified with the set (cf. Definition 4.3). By the invariance of 
vo\ k A(B) under rotation of p G S lfe , we get 



/ w\ k A(B) m - k+1 dW k = O k [ vo\ k A(B) m - k+1 d(S k - 1 ^ +i 
Jw k JM k 



Note that, up to a scaling factor, the right-hand side above is the coefficient 
C(m, k) introduced in Definition 4.3. Using (O k /O m )G k +i,m+i = Gk,m we obtain 



vo\K k (e) 



Qn 



C{m,k) f t m - k r(t) km - 2 X m (t) n - k - 1 dt. (21) 
■/ o 



This finishes the proof. □ 



Proof of Lemma 4-6- For given p G £ fc , let S(p ± ) denote the (k— l)-subsphere of 
S k perpendicular to p. At a point (B,p) G Wfc we have <iWfc = dS(p- L ) k+1 A cLS^ 
and we have dS'(p- L ) fe+1 = ^(p- 1 ) A • • • A dSip 1 ) at the point B = (b 1: . . . , b k+l ). 
The Jacobian J(.A) we are looking for is hence determined by 

d(S k ) k+1 (A) = J {A) dSip^f^iB) A dS k {p) A dt. 

Choose an orthonormal moving frame e±, . . . ,e k , e k+ \ with p{A) = e k+ \. For 
1 < i < k define the one-forms \i{ := —{ei,dp) (compare §3.2). Then the volume 
form of S k at p is given by dS k {p) = fi± A ■ ■ • A [i k . 

Differentiating te = Ap we get e dt = dAp + A dp. By multiplying both sides 
with A" 1 and using Formula (14) we obtain 

p(dt/t) -dp = A~ x dAp. 

Let Q denote the (k + 1) x (k + 1) matrix having the as rows. With respect to 
this basis, the above equation takes the form 

(m,...,H k , dt/t) T = Q(p (dt/t) -dp) = Q A' 1 dAp. (22) 
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Wedging the entries on both sides yields 

dS k (p) Adt = tvo\(AY l {p,dai) A ••• A (p,da k+1 ). (23) 

The volume form of S(p ± ) at bi is given by 

dSip 1 ) = (ei, dbi) A • • • A (e fc _i, dbi). 

In order to compare dS(p- L ) k+1 AdS k Adt with d(S k ) k+l we use a different moving 
frame. Fix an i, 1 < i < k + 1, and choose the moving frame as above and 
additionally with = 6j. Consider the modified frame ei, . . . , e^+i that arises 
after rotating bi to ctj and leaving the orthogonal complement of span(aj, fixed, 
i.e., e.j = ej for 1 < j < k — 1, := a^, and e^+i = — tfej + rp (cf. Figure 4). 




Figure 4: The frame (e^) and its modification (e,). 



This implies (e& + i, dai) = r(p, dai) — t(bi, dai) = (l/r)(p, dai), where we have 
used that bi = (a^ — tp)/r for the last equality. Hence the volume form of S k at 
di equals 

dS k (cii) = (l/r)(e l ,da i ) A • • • A (e fe _i,doj) A (p,da,i). 

If we wedge the (ei,dai) A ■ ■ ■ A (ek-\,dai) to both sides of Equation (23), we 
obtain on the right hand side 

k-\-l 

' f\ (eudai) A • • • A (e^dai) A <p, da t ) = -i- r fe+1 d(S fc ) fc+1 (A). 



vol(A) / = \ N x ' " */ xr, vQ|(A) 
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On the left-hand side we obtain, using (ej, den) = r ( e j> dbi) + t{ e j,dp) and taking 
into account that (ej,dp) A dS k (p) = since dS k (p) = f\^ =1 (ej , dp) , 

fc+ifc-i fc+ifc-i 

/\ /\ (ej,dai) A dS k A dt = r ( fc " 1 )( fe + 1 ) j\ /\ ( ej , d&i) A dS fe A dt 

1=1 J = l j=l jr' = l 

= r( fc - 1 K fc+1 >dS(p 1 ) fc+1 A A dt. 
This implies that J {A) = t' 1 r( fc - 2 )( fc+1 )vol(^) as claimed. □ 



4.2 The infeasible case 

Recall that T n ,m = (S m ) n \J-'n,m denotes the set of infeasible instances. We define 
for / C [n] 

Zi,m( £ ) '■= [ A g Zn,m \ ^(A) > e^ 1 and A has a SIC with blocking set /}. 

We note that by symmetry, the volume of Z^ m (e) only depends on the cardinality 
of/. 

Lemma 4.7 l^ m (e) has measure zero if\I\ > m + 1. 

Proof. If A 6 2^ m (e), then {aj | i £ /} is contained in the boundary of a SIC of 
A with blocking set /. Hence the affine hull of (cij)i g / has dimension less than m. 
However, if |/| > m + 1, the latter dimension is almost surely m + 1. □ 

It is known [5, Prop. 4.2] that in the infeasible case, blocking sets have at 
least m + 1 elements. This fact, together with Lemma 4.7, implies that 

""^-(rn + l) volJ K%)- ( 24 ) 

As with Fn^m, for ease of notation, we write T™ m (e) := zj^ 1 '^). 

The inequality in Theorem 1.3 for the infeasible case follows immediately 
from (24) and the following proposition. 

Proposition 4.8 We have for e € (0, 1], 



Qn 



< C(m,m) f (I - t 2 )^ {l - \ m {t)) n ' m ~ l dt. 
Jo 
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Proof. Consider the projection 

V: 2£ m (e) - (S™)"* 1 , A' = (a 1; . . . ,a n ) ^ A = (a 1; . . . ,a m+1 ). 

To investigate the image and the fibers of V> assume A' 6 X™ m (e). Then there 
exists p G S 1 ™ and q G (tt/2, 7r] such that cap(p, a) is a SIC of A' G 2"™ m (e) 
with blocking set [m + 1]. Then we have that (oj, p) = t for all i G [m + 1] and 
(di,p) > t for all i 6 [n] \ [m + 1], where t := cos a G [—1, 0). Lemma 3.2 implies 
that tp G convjai, . . . , a m+ i}. In turn, Lemma 3.1 implies that cap(— p, ir — a) is 
a SIC of A with blocking set [m + 1] and we obtain that A G X"™ +1 m( £ )- These 
reasonings show that the image of ij) is contained in F^+i m( £ )- 

Suppose now that a±, . . . , a m+ \ are linearly independent. Then it follows from 
Lemma 3.3 that the vector p is uniquely determined by A. This implies that the 
fiber of A under tp is contained in {A} x cap(p, a) n_Tra_1 . We conclude that for 
almost all A G T^ +X m 

^From these observations we obtain by Fubini 

voIX™ (e) /" i,i 



/-/->n— m— 1 — / 



In the proof of Proposition 4.5 we derived from the integral representation (19) 
P- the Forr 

vol2£ m (e) 



f° r V °'C-'^- ( 'i £ ' > the Formula (21). In exactly the same way we can show that 



< C(m, m) f (1 - t 2 )^ (1 - A m (t)) n ~ m_1 dt, 
Jo 



Qn 

which proves the assertion. □ 

Remark 4.9 (i) It may be of interest to compare, in the case m = 1, the 
upper bound for p(n, I, a) which follows from our results with the exact 
expression (1) for this quantity shown by Stevens. Recall, the latter is 

**>.•>- »(-r , -G)(-^r-^ i ©0-^)""' 

where k = [^\. For a G [0,7r/2], Proposition 4.8 and Proposition 1.5 yield 

p(n,l,a) = Prob{^GX;,i} + Prob{A0X- nil and^(^) > ^-1 

[ cos(a) J 

(\ /* COS Ct 

l)c{i,i) J o (i-t 2 r 1/2 (i-Ai(t)r 2 dt 



< 
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We use now that C(l, 1) = §, as shown in Example 4.4. Then 




arccos t 



) 



n-2 



p(n,l,a) = 



dt 



IT 



Oi/'K 




That is, we get Stevens' first term. 

(ii) It may also be of interest to compare, for the case m = 2, our upper bound 
for p(n, 2, a) with the upper bound in (2) obtained by Gilbert [10]. Recall, 
the latter gives 



where A denotes the fraction of the surface of the sphere covered by the cap 
of radius a. It is easy to see that our bound implies 



The first term in this sum is negligible for large n. The second term com- 
pares with Gilbert's for moderately large caps but it becomes considerably 
worse for small caps. 

5 On the values of the coefficients C(m,k) 

In this section we provide estimates for the numbers C(m,k). In §5.1 we derive 
upper and lower bounds for them which are elementary functions in m and k. 
In the case m = k the upper bound is actually an equality, yielding an exact 
expression for C(m, m). Then, in §5.2 we use these bounds to prove Theorems 1.2 
and 1.4. Finally, in §5.3 we briefly describe how to derive an exact expression for 
C(m,m — 1) and how, for any given m, one may obtain the values of the C(m, k), 
k = 1, . . . , m, by solving an m x m linear system. 

5.1 Bounding the coefficients C(m, k) 

Our first result provides bounds for C(m, k) in terms of volumes of spheres. 



p(n, 2, a) < -n(n - 1)A(1 - A) 



n-1 



p(n, 2, a) < — (n 2 - n + 2) + -n(n - 2)(1 - A) 



n-1 
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Lemma 5.1 We have for 1 < k < m 

k + l Ok-lOm-k , n( , , , (fc + l) m ~ fc+1 Ofc-lOm-fe 

— ; ^ < C(m,k) < : — . 

with equalities if k = m. In particular, C(m, m) = °q~ 1 ■ 
Proof. Recall the Definition 4.3 of the C(m, k) 



C ^ fc ) = —im — s k , m / Hi A) m - k+1 ds. 

U m J M k 



We set S := (S k ~ r ) k+1 and denote by U the open dense subset of S consisting of 
all B = (&!,... , bk+i) such that every k of these vectors are linearly independent. 
By Definition 4.1, M k is contained in U. 

Set A(£?) = conv{6i, . . . , b k +i} and Aj(5) = conv(0, 6i, . . . . . . , fyfc+i) (where 

6j means that 6j is omitted). We define for B £ U 

k+1 



mvol(.B) := J^volAi(S). 



i=i 



For S 6 Mfc we clearly have mvol(i?) = volA(S), but in general this is not the 
case. 

The essential observation is now the following: 

{vo\A) m ~ k+1 dS = \ [ mvo\ m - k+1 dS. (25) 

M k 2 Ju 

In order to show this, note that for B G U there exists a unique \i G R k+1 
with = 1, fJ>\---fJ>k 7^ 0, and such that X^=i Mi^j = 0- This allows to 

define the map (ft: U — > {—1, l} fe , £> i-> (sgn(^i), . . . ,sgn(//fc)). Note that = 
</> _1 (l, . . . , 1). Moreover, each a € {—1, l} fc defines an isometry 

M fc -> _1 (cr), B aB := (o"i6i, . . . ,a k b k ,b k+1 ). 

It follows that mvol(B) = mvol(cjS) since changing the signs of rows does not 
alter the absolute values of determinants. This implies 

/ mvor- fe+1 d,S= V / mvor- fe+1 d,S = 2 k [ mvo\ m - k+1 dS, 

Ju CT e{-i,i} fc k 

which proves the claimed Equation (25). 
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Recall now the norm inequalities 



(x{ + ■ ■ ■ + x e p ) < (xi + ■ ■ • + x p ) e < p l ~\x\ + ■ ■ ■ + zj) for xi > 0,4 > 1 (26) 

where the last follows from the convexity of the function R — > R, y ^ y^ . For the 
upper bound in the statement we now estimate the right-hand side of Equation 
(25) using the last inequality above (with p = k + 1 and i = m — k + 1). We 
obtain 

/ mvo\ m - k+1 dS < (k + l) m ~ fe / (volAi) m - fc+1 ^ 
= (jfe + i) m - fe +! / (volA fe+1 ) m - fc+1 d5 

where B 6 K hfe denotes the matrix with rows &i, . . . , Since the integrand on 
the right does not depend on bk+i, we can integrate over B G (S' fc_1 ) fc and pull 
out a factor Ofc_i obtaining 

m— fc+1 r 

k-l\k 



We plug in here the formula of the moments from Lemma 3.5. Putting everything 
together, and using Gk,m = {Pm-k/Om)Gk,m+i-, the claimed upper bound on 
C(m, k) follows. The lower bound is obtained by doing the same reasoning but 
now using the left-hand side inequality in (26). 

In the case k = m upper and lower bounds coincide and we get equalities for 
C(m,m). □ 



Remark 5.2 In the case k = 1 the upper bound in Lemma 5.1 coincides with 
the value for C(m, 1) shown in Example 4.4. 

For the proofs of Theorems 1.2 and 1.4 we need a more explicit expression for 
the bounds on the C(m, k). We devote the rest of this section to deriving such 
expressions. 
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Lemma 5.3 For 1 < k < m we have 



Cfe-lCm-fc < /tT 3/4 //m 



Cm V 2 \l \ k ^ 

In the cases k = 1 or k = m one has the sharper bound 2( ^"~ 1 < Jm. 

The proof uses bounds on Gamma functions, which we derive next. 
Lemma 5.4 For all r > 1, 

Proof. The double factorials /c!! are defined as follows. For k even, k\\ := 
k(k — 2)(fc — 4) ■ ■ ■ 2, and for fc odd, fc!! := k{k - 2)(k - 4) ■ ■ ■ 3 ■ 1. Also, by 
convention, 0!! = 1. By the functional equation T(x + 1) = xT(x) of the Gamma 
function, it easily follows that for r G N, r > 1, 

r + l\ J y^(r - 1)!! 2 _I1 2 i if r is even, 

2 / \ (r- l)!!2- r 2 i if r is odd. 

We estimate now double factorials in terms of factorials. If r > 2 is even, 

(r!!) 2 = rr(r-2)(r-2)---4-4-2-2 

= r ( r -l)-^-(r-2)(r-3)^|---4-3^2-2 

= r \ r r ~ 2 4 2 
r — 1 r — 3 3 



r r r-2r-2 44 

2-2. 



r — lr — lr — 3r — 3 3 3 
We use that < for £ > 2 to deduce from this 

rWr + 1 < (r!!) 2 < r\V2r for r > 2 even. (28) 
Similarly, for r > 1 odd, one shows that 

(r!!) 2 =r! 



r r — 2 r — 2 5533 



1 r-1 r-3 r-3 442 2' 



which implies 



v . 



^i— < (r!!) 2 < rlv 7 ^ for r > 1 odd. (29) 

By applying the bounds (28) and (29) to (27) and noting that 2 1 / 4 < ^/f , the 
claim follows. □ 
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Proof of Lemma 5.3. Assume that 2 < k < m. Then, using Lemma 5.4, we 
deduce that 



o m r(|)r(^|±i) 

< 



(m — 1)! / m \ 1 ^ 4 



2 V (fc-2)!(m-ife-l)! 



7r f Tn\ (m — k)k(k — 1) / m 



2y \k J \ m \(k—l)(m — k) 



1/4 



The cases A; = 1 and k = m follow similarly from Lemma 5.4. □ 

An immediate consequence of Lemmas 5.3 and 5.4 are the following bounds 
on the coefficients C(m,k). 

Proposition 5.5 For 1 < k < m, 



In addition, for all m > 1, 



\ (m+ l)\fm 
C(m,m)<- !y . □ 



Remark 5.6 Using Lemmas 5.1 and 5.4 it is easy to obtain lower bounds for 
the C(m, k) similar to the upper bounds in Proposition 5.5. 

5.2 Proof of Theorems 1.2 and 1.4 

The following identity is repeatedly used in the proof. 

Lemma 5.7 We have J2n=k il) zn ~ k = ( l ~ z )~ k ~ l for At G N and z G (0, 1). 
Proof. Take the kth derivative on both sides of J2^=o z " = T=z - a 
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Proof of Theorem 1.2. By definition, we have N(m, a) > n iff cap(ai, a) U • • • U 
cap(a ra ,a) / S n . Hence 

oo oo 

E(iV(m, a)) = ^2 Prob(iV(m, a) > n) = y~]p(n, m, a). 

ra=0 n=0 

We assume that a < tt/2. Since p(n, m, a) = 1 for n < m + 1, we conclude 

oo 

E(iV(m,a)) = m + 1 + ^ p(n,m,a). (30) 

n=m+l 

Proposition 1.5 states that, for a G (0, |], and £ = cos(a) 



p(n, m, a) = 




where we have put 

Pn,m( £ ) ■= P™b{A G J„, m and <g(A) > e" 1 }- 
We first estimate 

E 2l ~" E (" ~ k ) 

n=m+l fc=0 v 7 

as follows (take r = n — 1) 

r = E E (I.) Q) Q) E I G 

fc=0r=m \ / \ / fe=0 v 7 r =fc V 7 V 7 fc=0 

Applying Lemma 5.7 to the last expression we obtain 

T < y ( i 2 fc+1 - 2 + — E_ < 2m + 1. 

fe=o v 7 

We now estimate T* := E^^l m+1 P n ,m{£) using Theorem 1.3, which tells us 
that 

Pn, m (e) < (j % + ^C(m,m)J\l-t 2 ) m ^(l-\ m (t)) n - m - 1 dt. 
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Hence, using Lemma 5.7 again, 

^ \ re 2 „ 

\n-m~l dt 



T* = C(m,m) J] L + J / (l-t 2 )^(l-X m (t)y 

n=m+l ^ /JO 
re 00 



n=m+l 



m + 1 



C{m,m) ( \ m (t)- rn - 2 dt < C{m,m) 
Jo 



.ra+2 ' 



Plugging in the estimate for C(m, m) from Proposition 5.5, we obtain the claimed 
bound for E(iV(m,a)) < m + 1 + T + T*. □ 

We now turn to Theorem 1.4 on the expected value of lnff^A). In Theo- 
rem 1.3 we derived tail estimates on the probability distribution of the GCC 
condition number. For the sake of clarity, we include the following simple obser- 
vation showing how to use these tail estimates to bound the expected value of 
the logarithm of the condition number. 

Proposition 5.8 Let Z be a random variable, almost surely greater or equal 
than 1, satisfying, for some K,to > 0, that Prob{Z > t} < Kt^ 1 for all t > to- 
Then E(lnZ) < lnt + f . 

Proof. We have ProbjlnZ > s} < Ke~ s for all s > Into- Therefore, 

f'OO roc 

E(lnZ)= / Prob{lnZ > s}ds < lnt + / Ke~ s dt = lnt + — , 

Jo ./Into *0 

as claimed. □ 
We next proceed to prove Theorem 1.4. To simplify notation we put 

Pn,m(e) ■= Prob{A G J n , m and V(A) > e' 1 }, 
Qn,m( £ ) ■■= Prob{A G ^ n , m and ^(A) > e- 1 }- 

Lemma 5.9 For any n > m > 1 and e G (0, 1] we have 

(i) If e- 1 > 13(m + 1) 3/2 then P„, m (e) < 2e(m + l) 3 / 2 e. 

(ii) Jfe" 1 > (m + l) 2 then Q n , m (e) < (m + 1) 7/4 e. 
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Proof, (i) Theorem 1.3 tells us that 

Pn,m(e) < ( n + J C(m,m) - t 2 )^(l - A m (i))"— 

Recall Formula (5) for the relative volume A m (t) of a cap of radius arccos(t) 
on S m . Recall also from Lemma 5.3 that a m := 2( ^ n ~ 1 < Jm. The first order 
derivative of X m (t) 

dA^t) 1 2 ^ 

(ft 2 v 7 

is increasing, hence A m is a convex function. Moreover, A m (0) = 1/2. This 
implies 2\ m (t) > 1 — a m t for all t G [0, 1], hence 1 — A m (i) < ^(1 + a m t). 
Bounding C(m,m) as in Proposition 5.5 we arrive at the estimate 

Pn, m (e) < 2(m + 1) ^ ( m " x ) (1 + v^r™- 1 e. (31) 

We now proceed dividing by cases. Suppose that e~ l > 13(m + l) 3 / 2 . 
Case 1: n < 13(m + 1). In this case > n^/m and hence, using (31), 

Pn,m(e) < 2(m + 1) Vm (1 + l/n)"-"^ 1 e < 2e(m + l)Vm e. 

Case 2: n > 13(m + 1). This implies In (e^^j < In (|), and it follows 
that 

n \ 1 , L , / en \ / 4 \ 



m + i7 £ (^TI)I' ! "' +l£ l^ £ (5j- <32) 

Since, in addition, e _1 > 13(m + l)yfm > 2yjm we get from (31) 

P n , m (e) < 2(m + ^ ( m " J (0 e < 2(m + 1) £• 

(ii) Theorem 1.3 implies that 



m / \ 

< E ( fc + ij C ( m ' fcje"'-^^-^-*- 1 ) 

< ^C(m,fc)e m - fe+1 2 fe+1 , 



rt— fc— 1 



(ft 



fc=l 
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the second line since \ m (t) < \. Using Proposition 5.5 we obtain 



< £ ^(m + l) 7 / 4 ^(^) ((m + l)e) m - fc 

k=l ^ ' 

< £^( m + l) 7 / 4 (l + (m + l)e) m 

Under the assumption e^ 1 > (m + l) 2 we have (m + l)e < and hence 

Q„, m (e) <£\^(m + l) 7 / 4 ^. □ 

Proo/ o/ Theorem 1.4- For > 13(m + l) 2 we have by Lemma 5.9 
Probj^) > e- 1 } = P n , m (e) + Q„, m (e) 

< (2e(m + 1) 3/2 + v 77 ^^ + l) 7/4 )e 

< 9.6(m + l) 2 e. 

An application of Proposition 5.8 with i^T = 9.6(m + l) 2 and to = 13(m + l) 2 
shows that 

E(lntf(A)) < 21n(m+ 1) + In 13 + 9.6/13 < 2 ln(m + 1) + 3.31. □ 
5.3 On calculating the C(m, k) 

We describe a method for calculating the C(m, k). For 1 < K m < n we define 
the following integrals 

J(n, m, fc) := 2™- 1 ^ ™ J * t m " fc (l - i 2 )**™"^*)"-*- 1 ^. 

By setting e = 1 in the first part of Theorem 1.3 we get from (4) that, for all 

n > m, 

m m—l , _. s 

Y J An^k)C{m,k) = Y J \~ k )• (33) 
k=l k=0 ^ ' 

By taking m different values of n (e.g., n = m + 1, . . . , 2m) one obtains a (square) 
system of linear equations in the C(m, k). Solving this system with Maple (sym- 
bolically for even m and numerically for odd m) we obtained Table 1. 

We can further use (33) to obtain expressions for C(m, k) for values of k other 
than 1 and m. We do so for k = m — 1. 
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k\m 


1 


2 


3 


4 


5 


6 


1 


■2 

7T 


2 


5.0930 


12 


27.1639 


60 


2 




3/4 


3.9317 


477/32 


49.5841 


78795/512 


3 






0.6366 


39/8 


25.1644 


897345/8192 


4 








15/32 


4.8525 


132225/4096 


5 










0.3183 


4335/1024 


6 












105/512 



Table 1: A few values for C(m, k). 



Proposition 5.10 For all m > 2, 

C(m,ro-1) = 2m _ 1 (l + O where a m = — . 

Sfcefc/i o/ Proo/. Put J(n, m, fc) := £ t m " fc (l - i 2 )5 fem - 1 A m (i)"- fe - 1 cii so that 
I(n,m,k) = 2 n - 1 ( fc ™ 1 ) J(n,m, fc). In the following we write N = n — m. One can 
prove that for fixed m the following asymptotic expansion holds for N — > oo 

on-m-i t / \ 11 m(m - 1) 1 

2" m V^m,™) = — — - 



2 n ~ m J(n,m,m- 1) 



a m , iV 
1 



a 3 iV 3 



a2 1 (iV + l)(iV + 2) + U 4 



It follows after a short calculation that the left-hand side of (33) has the following 
expansion 



C(m, m) 



—N rn + "'^ r- 1 + ( - m ( m ~ jj ] 7V m ^ 2 



(m + 1)! \ 

2"i— i ^ 

+ C(m, m — 1) 



■N 



m-2 



+ O (N m - 3 ) 



where 



ai(m) := ^ j' = -m(m + l), 

0<j<m 

a2(m) := ij = — (m + l)m(m — l)(m — 2)(3m + 2). 



0<i<j<m 
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Now we expand the right-hand side of (33) to obtain 
ml \ ml (m — 1)1 J 

/ *MI + + o<jv— ). 

\ m! (m — 1)! (m — 2)!/ 

By comparing the coefficients of N m (or those of _/V m_1 ) we obtain 

C(m,m) = 2m a m . 

By comparing the coefficients of N m ~ 2 we get, after a short calculation, 

O 2 _ ( 

C(m,m-1) = ^,2 2^-3 (^2 (ra ~ 1) ~ «2 (m) + max (m - 1) 

, m(m — 

+ m(m-l)+ V 

m— 1 / 

and simplifying this expression, the claimed result follows. □ 
Finding a closed form for all coefficients C(m, k) remains a challenging task. 
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